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Abstract 



In this article we study the geodesic motion of test particles and light in the five-dimensional (ro- 
tating) black string spacetime. If a compact dimension is added to the four-dimensional Schwarzschild 
or Kerr spacetime, the new five-dimensional metric describes a (rotating) black string. The geodesies 
in the Schwarzschild and Kerr spacetime have been studied in great detail, however, when a compact 
dimension is added new behaviour occurs. We present the analytical solutions of the geodesic equa- 
' , tions and discuss the possible orbits. The motion in the ordinary four-dimensional Schwarzschild and 

^jH- Kerr spacetime is compared to the motion in the (rotating) black string spacetime. 

U 

i^B- 1 Introduction 

In 1916 Schwarzschild introduced the first exact solution of the vacuum Einstein equations in four dimen- 
sions pp. This unique solution describes a spherically symmetric black hole. About fifty years later, in 
• 1963, the solution of a rotating black hole was discovered by Kerr [2J. 

If a compact dimension is added to the Schwarzschild metric, it describes a five-dimensional black 
^sO . string. Analogously the Kerr spacetime turns into the spacetime of a rotating black string if a compact 

dimension is added. 

The idea of compact dimensions came up in the 1920s after Kaluza suggested to add a fifth dimension 
to general relativity [3]. This was one of the first attempts to unify gravity and the electromagnetic forces. 
In 1926 Klein proposed a physical interpretation pi]: The added dimension is compact, i.e., it is curled up 
in itself and has a certain length, which is too small to be observed. But still the theory of Kaluza and 
Klein could not explain the weakness of gravity in comparison to the electromagnetic forces. 

In the 60s and 70s the string theory was developed. Here the existence of higher dimensions is essential 
for the internal consistency of the theory, mostly they are also assumed to be compact. 

The equations of motion for test particles in the Schwarzschild spacetime were solved in 1931 by 
Hagihara [5] in terms of the elliptic Weierstrafi p-, a- and ^-functions. In the Taub-NUT [5J, Reissner- 
Nordstrom [7[ and Myers-Perry [5] spacetime the equations of motion were also solved in terms of the 
elliptic Weierstrafi functions. 

The integration of geodesies was advanced in the papers of Hackmann and Lammerzahl [51 [TU] . They 
integrated the geodesies in the four-dimensional Schwarzschild-de Sitter space-time analytically in terms 
of the hyperelliptic 9- and cr-functions. The mathematical problem is based on the Jacobi inversion prob- 
lem which can be solved if restriced to the 0-divisor. This method was also applied to find the analytical 
solution of the equations of motion in higher dimensional Schwarzschild, Schwarzschild-(anti)de Sitter, 
Reissner-Nordstrom and Reissner-Nordstrom -(anti) de Sitter spactime [11] as well as in Kerr-(anti) de 
Sitter spacetimes [T2J and in higher dimensional Myers-Perry spacetime [13] . Recently, the geodesies equa- 
tions in special cases were solved analytically in the singly spinning black ring spacetime |14| and in the 
(charged) doubly spinning black ring spacetime [TS] . 

The test particle motion in different spacetimes containing black strings was studied in [TBI H3 [HI UHl 
[20]. In [2TJ and [22] the analytical solutions of the equations of motion in the Schwarzschild and Kerr 
spacetime pierced by black string were presented. Furthermore, the geodesic motion of test particles in 



field theoretical cosmic string spacetimes were investigated, namely Abelian-Higgs strings |23j . two inter- 
acting Abelian-Higgs strings (51] and cosmic superstrings [25] . 

In this paper we will consider static and rotating black strings, which can be obtained by adding an 
extra compact dimension to the Schwarzschild and Kerr metric. The geodesic motion in the ordinary 
Schwarzschild spacetime and the Kerr spacetime has already been analyzed in detail, however, when 
a compact dimension is added new behaviour occurs which is studied in this article. We present the 
analytic solutions of the geodesic equations in terms of the elliptic Weierstrafi functions and discuss the 
corresponding orbits. In the first part the non-rotating black string is studied and compared to the 
Schwarzschild black hole. In the second part the rotating black string is analysed and compared to the 
Kerr black hole. 



2 Black string spacetime 

In this section we will discuss the geodesies in the static black string spacetime and present analytical 
solutions of the equations of motion. 

2.1 The geodesic equations 

The well known Schwarzschild metric is the unique static spherically symmetric solution to the vacuum 
Einstein equations in four dimensions. If an extra compact spatial dimension w is added, the metric takes 
the form 

2M\ .of 1 



da' = - ^1 - — J dr + y- 2mJ dr 2 + r*(6& + sin' M<p*) + dw 2 , (1) 

where M is proportional to the mass of the black hole. This solution describes a neutral uniform black 
string. As usual the singularity is located at r — 0. The horizon is located at r = 2M and covers the 
extra dimension. The gravitational constant and the velocity of light are set to 1. 
The Hamilton- Jacobi equation 

!fr + 2 s " d^d^ = ° ^ 

has a solution of the form 

S=-ST-Et + Ltp + Jw + S r (r). (3) 

Here r is an affine parameter along the geodesic. The parameter 5 is equal to 1 for particles and equal to 
for light. E, L and J are the conserved momenta; E is the energy, L denotes the angular momentum 
and J is a new constant of motion according to the compact dimension w. We set ■& — 5 since the orbits 
lie in a plane if plotted in Cartesian x-y-z coordinates (x = r cost/? sin i?, y = r sintpsinz?, z — rcosi?) due 
to the spherical symmetry of the original Schwarzschild metric. 

For convenience, we introduce dimensionless quantities (rs = 2M) 

r=— , w = — , t = — , t = — , L = — . (4) 

rg rs rs rs rg 

In these dimensionless quantities the horizon is at fn = 1- 

The Hamilton- Jacobi equation ([2]) separates and yields a differential equation for each coordinate 



,~\ 2 

dr 



d~, 

dtp 

d") 
dw 
d-, 

di Ef 3 



dj f — 1 

with the polynomial 



R (5) 
L (6) 
Jr 2 (7) 

(8) 



R = E 2 ? 4 -(5 + J 2 )(P - p) - L 2 {f 2 - f) . (9) 



We also used the Mino time 7 as f 2 d^j — df [26 . The new constant of motion J due to the compact 
dimension is only present in the f -equation (JSJ and of course in the iD-equation ([7]) . The other equations 
are the same as in the original Schwarzschild spacetime without the compact dimension. 

2.2 Classification of the geodesies 

Before discussing the f -motion in detail, we introduce a list of all possible orbits: 

1. Terminating orbit (TO) with ranges f £ [0, 00) or r £ [0,ri] with r\ > ?h- The TOs end in the 
singularity at f = 0. 

2. Escape orbit (EO) with range r € [r*i,oo) with r\ > fjj. 

3. Bound orbit (BO) with range r E [r\,r2\ with t\ < ri and t\,t-z > fn- 

The right hand side of the differential equation ((SJ) has the form R = ^ i=1 onf 1 with the coefficients 



a 4 = E 2 — 5 — J 2 (10) 

a 3 = 6 + J 2 (11) 

a 2 = -L 2 (12) 

ai = L 2 . (13) 



R > is required in order to obtain real values for f. The regions for which R > are bounded by 
the zeros of R. The number of zeros can be determined with the help of the rule of Descartes. One can 
distinguish four cases: 

1. If E 2 > 6 + J 2 and L 2 > 0, then two or no zeros are possible. 

2. If E 2 > 5 + J 2 and L — 0, no zeros are possible. 

3. If E 2 < S + J 2 and L 2 > 0, then one or three zeros are possible. 

4. If E 2 < 5 + J 2 and L — 0, one zero is possible. 

Since E is real and therefore E 2 > 0, S + J 2 > is required. So if J = (where the equations of motion 
reduce to the original Schwarzschild case without the compact dimension) the cases 3. and 4. are not 
possible for light (S = 0). But if J 7^ there are 3 possible zeros of R for 5 — 0, that means in the 
Schwarzschild spacetime with an extra compact dimension we have the possibility of stable bound orbits 
for light, which are not possible in the original Schwarzschild spacetime. 

We define an effective potential V from the equation §5§ by 

(^y = f\E 2 -V), (14) 

thus 

In the limit f —} 00 the effective potential V converges to 5 + J 2 . If f —> 0, we have V — > —00. (5^) = 
determines the turning points of an orbit. Depending on the number of turning points and the shape of 
the effective potential four different types of orbits are possible: 

• Type A: no zeros. Only TOs are possible. 

• Type B: one zero. Only TOs are possible. The turning point can coincide with the event horizon. 

• Type C: two zeros. TOs and EOs are possible. If the energy E coincides with the maximum of the 
effective potential then there is a single zero, which corresponds to an unstable circular bound orbit. 

• Type D: three zeros. TOs and BOs are possible. 

If the energy E coincides with the minimum of the effective potential then there are two zeros. TOs 
and circular BOs are possible. 

Table [1] shows an overview of the types of orbits. Some examples for energies corresponding to certain 
orbits in the effective potential can be seen in figure [TJ 



type 


zeros range of f 


orbit 


A 


o ' 1 


' TO 


B 
B 


1 ' 1 • 

' * 


* TO 

TO 


C 


2 ' 1 • - 


TO, EO 


D 


3 ' 1 • 


• ' TO, BO 



Table 1: Types of orbits of light and particles in the black string spacctimc. The thick lines represent the range of the orbits. The 
turning points arc shown by thick dots. The horizon is indicated by a vertical double line. 
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(a) 8 — 1. L — >/2 and J — 1.5: Examples of the (b) 8 — 0, L — \/5 and J — 1: Examples of the 
orbit types A, B and Bo- orbit types C and D. In the case of type D, 

here a stable bound orbit for light is possible. 

Figure 1: Examples of the effective potential in the black string spacctime. The grey area is a forbidden zone, where no motion is 
possible. Horizontal red dashed lines represent energies and red points mark the turning points. The horizon is marked 
by a vertical black dashed line. 



2.3 Solution of the geodesic equations 

In this section we present the analytical solution of the equations of motion ([5])-([H])- 



2.3.1 The f -equation 

The polynomial R = Y^t=i can be reduced to cubic order by the substitution f — ±-: R' = Y^=o biX 1 . 
A further substitution x = ^ (Ay — y) transforms R' into the Weierstrafi form so that equation ([5]) turns 
into 



dry 



4y J - 522; - 53 = P 3 (y) , (16) 



where 

_ h\ hb 3 _ bib 2 b 3 b b 2 3 b\ 

92 ~ 12 " ~T ' 93 ~ ~8 16" ~ 216 ■ 
The differential equation (fT6)) is of elliptic type and is solved by the WeierstraB p-function [27] 

2/(7) = Pil ~ 1^92,93) , (18) 
where 7/ = 7; n + f°° . dy with y; n = ±J£- + If. Then the solution of JS]) acquires the form 

J 2/in ^4y 3 -g 2 y-g 3 4r ' n 



b 



± 1" ■ ( 19 ) 

4p(7-7i' n ;52,53) - f 



2.3.2 The (^-equation 

Equation §§§ has the trivial solution 



99(7) = £(7 - Tin) + <£in 



(20) 



2.3.3 The ^-equation 

Using equation (J5j> and the subsitution f = ± - ba b2 the iD-equation (0 becomes 

d ^= Jhl - 2 7 a = . (2D 

The right hand side is an elliptic differential of the third kind. After the substitution y — p(v) with 
p'(v) = ^/4p 3 (v) — Q2p{v) — 33 and integration of equation (|2Tj) one gets 

«(7) - = ^ / — dv> , (22) 

lb J Vin (p(v') - p{u)) 

where v = v(j) — 7 — 7i„,i>i n = u(7in) an d p(u) = f§ = — fj- ^he s °l ut i° n °f this integral is given in 
terms of the elliptic p-, er- and ^-function [51(7]: 



. J6 2 . p"(u) /„,., w . c(u — u) , cr(u + u) 
w ^) = -^ilbhk (2CW(^-«m)+ln-A -L-In 



16 (p'(u)) 3 V " (T (u in - li) ^(Uin+u), 



2fj(tt)(w - Win) + 2(C(U) - C(Uin)) + 



(23) 



16 (p'(u)) 2 V ' " p{v)-p{u) p(v in )-p{u) 

2.3.4 The t-equation 

Using equation ([5]) the f-equation (j8]) becomes 

Ef 3 df . , 

dt = ■ . (24) 

Now we substitute r — ± — 2a__ anc [ apply a partial fractions decomposition. Then we get 

4y 3- 

T 2 T" 2 77 1 7" 2 771 r 2 

where pi = — j^, P2 = Ai = and = — The right hand side consists of elliptic differentials 
of the third kind. After the substitution y — p[v) with p'(v) — ^/4p 3 (w) — g2p{v) — <?3 and integration of 
equation (|23|) one gets 

f(7) " im = L (S pM-W) + (R^wf) dv ' ' (26) 

where w = ^(7) = 7 — 7; n ,^in = ^(7in) and pi = p(vi). The solution of this integral is given in terms of 
the elliptic p-, a- and (- function OH]: 

*(7) = >, "TT^ 2 C0w - w in) + In -7 r - In — — — 

~^P{Vi)\ cr(v in -Vi) u{v- m + Vi)J 

(p'(wi)) 3 \ cr(u in -i;i) u{v- m + vi)J 

Kl fo U XI \ 1 \ u P» P>in) \ . r 

2p(wi)(u - ^n) + 2(C(w) - C(vin)) H t\ ? — T 1 ^ 7 — \ + *in • 



(p'ivi)) 2 V p{v)-p{vi) p(w in )-p(ui) 



2.4 The orbits 



With these analytical results we have found the complete set of orbits for light and test particles in the 
black string spacetime. Depending on the parameters 5, L, J and E, BOs, EOs, and TOs are possible. 
BOs can be seen in figure [5] and [3] An EO is depicted in figure 0] and figure [5] shows a TO. In Cartesian 
x-y-z coordinates the orbits lie in a plane. However, this is not the case if the orbit is plotted in x-y-w 
coordinates. 

In the black string spacetime BOs for light are possible (see figure [5]). Such orbits do not exist in the 
ordinary four-dimensional Schwarzschild spacetime. 




(a) x-y-z-\Aot (b) x-y-w-\Aot 



Figure 2: 8 = 0, L = s/E, J = 1 and E = 0.973: 

Bound orbit for light in the black string spacetime. This orbit is not possible in the ordinary four-dimensional Schwarzschild 
spacetime. The sphere or cylinder is the horizon. 




3 Rotating black string spacetime 

In this section we will discuss the geodesies in the rotating black string spacetime and present analytical 
solutions of the equations of motion. The rotating black string metric is derived from the Kerr metric by 
adding an extra compact dimension. A detailed analysis of the geodesies in the original Kerr spacetime 
without a compact dimension can be found in e.g. |28j . 



(a) x-y-z-plot 



(b) x-y-w-plot 




Figure 5: 5 = 1, L = 2, J = 0.1 and _B = 1.0025: 

Terminating orbit for particles in the black string spacetimc. The sphere or cylinder is the horizon. 



3.1 The geodesic equations 

The Kerr metric is the rotating axially symmetric solution of the Einstein equations in four dimensions. 
If an extra compact spatial dimension w is added, the metric takes the form 

ds 2 = -%(dt - a sin 2 tfd^) 2 + ^V 2 + a 2 )d^ - adt) 2 + ^dr 2 + pW + du, 2 (28) 
p z p z ' A r 

where A r = r 2 — 2Mr + a 2 and p 2 = r 2 +a 2 cos 2 i?. This solution describes a rotating uniform black string. 
M is proportional to the mass of the black string and a is proportional to the angular momentum. The 
gravitational constant and the velocity of light are set to 1. The (ring) singularity is located at p 2 = 0, 
i.e. at r = and ■& = \. Hence geodesies with r — and i? ^ \ do not end in the singularity. So in the 
rotating black string spacetime negative r- values are allowed as in the Kerr spacetime |28) . There are two 
horizons defined by A r = and for a 2 < | they are given by 



r±=M± \/M 2 - a 2 . 



(29) 



Note that the metric is given in Boyer-Lindquist coordinates r, z9, tp. Cartesian coordinates can be obtained 
by the transformation 



= v (r 2 + a 2 ) sin $ cos ^ 



The Hamilton- Jacobi equation 



has a solution of the form 



y = \J (r 2 + a 2 ) sin ■& sin tp (30) 
z = r cos $ . 



1 



5= -Sr-Et + Ltp + Jw + S r (r) + S#(<&) . (32) 

Here r is an affine parameter along the geodesic. The parameter 5 is equal to 1 for particles and equal to 
for light. E, L and J are the conserved momenta; E is the energy, L denotes the angular momentum 
and J is a new constant of motion according to the compact dimension w. 
For convenience, we introduce dimensionless quantities (rg = 2M) 

r = — , w = — , t = — , t = — , L = — , a = — , K = — . (33) 
r s rs rs rg rg rg rg 

K is the famous Carter constant resulting from the separation of the Hamilton- Jacobi equation, see [29 . 



In these dimensionless quantities the horizons are located at f± = | ± y j — a 2 . 

The Hamilton- Jacobi equation (J5|) separates and yields a differential equation for each coordinate 

df x 2 
ch, 

d7 

^ = JL[(f 2 +c?)E-al] l —(hEsin 2 d - L) (36) 

d7 A r sin $ 

dw 

dj 

dt 



= R (34) 

f)' - e (35) 



Jp 2 (37) 
[(r 2 + a 2 )E - hL] - a(aE sin 2 - L) (38) 



d7 A 
with the polynomial R and the function 0: 

R = [(f 2 + h 2 )E - hL] 2 - K r {K +{S + J 2 )r 2 ) 

1 (39) 

6 = K- (S + J 2 )h 2 cos 2 2— (aEsm 2 d-i) 2 . 

sin ■& 

We also used the Mino time 7 as p 2 d7 = df [5^. The new constant of motion J due to the compact 
dimension is present in the f-equation (j34]) , the i9-equation (1351) and of course in the w-equation (j37]) . The 
other equations are the same as in the original Kerr spacetime without the compact dimension. To obtain 
the geodesic equations of the original Kerr spacetime, one can choose J = 0. 



3.2 Classification of the geodesies 

The function and the polynomial R define the properties of the orbits depending on the parameters of 
the metric and the test particle. In this section we will analyze the function O and the polynomial R to 
determine the possible types of orbits. A similar analysis of the Kerr spacetime can be found in [28], for 
the Kerr-(anti-) de Sitter spacetime see [12] . 

To obtain real values of 1? and f the conditions & > and R > have to be fullfilled, otherwise no 
motion is possible. From > we can immediately conclude that K > 0. If K — 0, the motion takes 
place in the equatorial plane [28] . 



Furthermore we see that a geodesic hits the singularity at p 2 = if K — (aE — L) 2 : A solution of 
p 2 = f 2 + a 2 cos 2 i? = is f = and simultaneously $ = ? . Since 

9 (-) = # - - Z) 2 > and 

V2^ V ; " (40) 

i?(o) = -a 2 [if - (as - l) 2 ] > o , 

it follows that K = (aE — L) 2 . Additionally an orbit lies in the equatorial plane (§ = ^) if K = (aE — L) 2 
(compare |12[ \2~E\ ). From (f4"U)) we also see the following: 



K > (aE — L). The geodesies cross i? = § but do not cross f = 0. 

K < (aE — L) 2 : The geodesies do not cross $ = \ but it is possible to cross f = 



Here the rotating black string shows the same properties as the Kerr black hole (see |28|). 
3.2.1 The ^-motion 

As in [T^] we substitute v — cos 2 d (with v e [0, 1]) in the function 0($): 

i? 

Q(y) =K-(5+ J 2 )~a 2 v - a 2 E 2 (l - v) + 2aEL - - — — . (41) 

Now we want to determine the number of real zeros of Q(v) in [0, 1]. The number of zeros only changes 
if a zero crosses or 1, or if a double zero occurs, v = is a zero of if 

Q(u = 0) = K- (aE-L) 2 = (42) 

and therefore 

L = aE±\[k . (43) 

Since v = 1 is a pole of Q(v) for L ^ 0, it is only possible that v = 1 is a zero of Q(is) if L = Q. 

Q(v = l,L = 0) = K- (5 + J 2 )a 2 = (44) 

yields that k = (5 + J 2 )a 2 . Note that due to the new constant of motion J (the momentum in the 
ill-direction) the value of K is different here in comparison to the Kerr spacetime. 

So L = aE ± \f~K and L = A K = (5 + J 2 )a 2 are border cases of the t?-motion. 

Next we consider the polynomial 

<d'(v) = a 2 (5 + J 2 — E 2 )v 2 + [2aE(aE — L) — (5 + J 2 )a 2 -K]v + K- (aE - L) 2 (45) 

to remove the pole of 0(v) at u = 1. Q and 9' are related by <d(i/) = j^Q'(v). Double zeros of @(i>) in 
[0,1) and hence of Q'(v) occur if 

9»=0 and ^M= . (46) 

From these conditions we obtain 



f E±(a 2 (8 + J 2 )-K)^E 2 -(5 + J 2 ) 

L 2a(5 + J 2 ) ■ (47) 

Therefore double zeros are only possible if E 2 > 5 + J 2 . 

With the help of equation (|43|) and (|47| parametric L-£' 2 -diagrams can be drawn. Figure |6] shows a 
typical example of a _L-_E 2 -diagram in the rotating black string spacetime, which resembles the parametric 
diagrams in the Kerr or Kerr-(anti-) de Sitter spacetime [12]. Here we recognize four regions. In region 
(a) and (d) we have Q(v) < for all v £ [0, 1] and therefore no geodesic motion is possible. In region (b) 
Q(v) has one real zero in [0, 1], here K > (aE — L) 2 and the orbit crosses d — ^. In region (d) Q(v) has 
two real zeros in [0, 1], here k < (aE — L) 2 and $ = \ is not crossed. However, crossing f = is possible 
in region (d) but not in (b). 

In the special case L = there is a single zero at vq = ^i^^^^j-iy^ (not necessarily in [0,1]). If 

additionally K = (6 + J 2 )a 2 then this zero is at v = 1 and if K = a 2 E 2 this zero is at v = 0. In the case 
K > (6 + J 2 )a 2 there is a zero v a in [0, 1] if K < ~a 2 E 2 and 8(^) > for all v e [0, 1] if K > a 2 E 2 . In the 
case K < (S + J 2 )a 2 there is a zero v a in [0, 1] if K > a 2 E 2 but Q(v) < for all v e [0, 1] if K < a 2 E 2 . 




Figure 6: 5 = 1, a = 0.45, J = 1, K = 2; 

Parametric Z-i? 2 -diagram of the tf-motion. In the grey regions (a) and (d) we have < for all v £ [0, 1] and hence 

no motion is possible. In region (b) O(f) has one real zero in [0, 1] and in region (d) 0(i/) has two real zeros in [0, 1]. 



3.2.2 The f-motion 

The zeros of the polynomial R are the turning points of orbits of light and test particles and therefore R 
determines the possible types of orbits. In contrast to the non-rotating black string, here we can consider 
negatives f- value too, since there is no singularity for r = and $ 7^ Orbits that cross f = are called 
crossover orbits [T2]. Before discussing the f-motion in detail, we introduce a list of all possible orbits: 

1. Transit orbit (TrO) with range r G (— oo,oo). This orbit is a crossover orbit. 

2. Escape orbit (EO) with range f € [r\, oo) with r\ > f+, or with range f £ (— oo, r\] with n < 0. 

3. Two-world escape orbit (TEO) with range [r%, oo) where < r\ < r~. The TEOs cross both horizons 
twice and emerge into another universe. 

4. Crossover two-world escape orbit (CTEO) with range [ri,oo) where n < 0. The CTEOs cross both 
horizons twice and emerge into another universe, f = is crossed once. 

5. Bound orbit (BO) with range f £ [ri,T2] with < n < r-z and 

(a) either r± , r% > r + or 

(b) n , r 2 < r_ . 

6. Many-world bound orbit (MBO) with range f € [^i,^] where < r\ < r_ and r2 > r + . The MBOs 
cross both horizons several times. Each time both horizons are traversed twice the test particles 
emerge into another universe. 

7. Terminating orbit (TO) with ranges f € [0, oo) or f G [0, r±\ with 

(a) either r% > f + or 

(b) < n < f_. 

The TOs end in the singularity at f = and i? = ^. AT = (a_E — Z) 2 is required. 

The possible orbit types depend on the number of real zeros of the polynomial R. The number of zeros 
can only change if double zeros occur, i.e. if 

R(r)=0 and ^^ = 0. (48) 
dr 

With the help of equation (|4"5|) parametric L-i? 2 -diagrams can be drawn, where five regions with a different 
number of real zeros of R appear (see figure [7]). Also the ^-equation has to be taken into account. Both 
parametric L-i? 2 -diagrams for the ^-motion and the f-motion have to be combined in order to detect all 
possible behaviour of light and test particles in the rotating black string spacetime (see figure [SJ . 




Figure 7: 5 = 0, a = 0A5, J = 1.3, K = 2: 

Parametric L-E 2 -diagram of the f-motion. In region (I) R has no real zeros. In the regions (II) and (V) R has two zeros. 
Four zeros arc possible in the regions (III) and (IV). Here we see that in the rotating black string spacctimc all five regions 
arc present for (5 — 0, whereas in the Kerr spacctimc region (III) and (V) do not exist if 5 — 0. 



L 




(a) S = 1, a = 0.45, J = 1, K = 2 (b) S = 1, a = 0.45, J = 1, K = 2: (c) 5 = 1, a = 0.45, J = 0.1, K = 0.2: 

Different view of figure (a), here region For small K a second part of region 
(Id) is visible. (IVb) appears. 



Figure 8: Combined parametric L-i? 2 -diagram of the i?-motion and the f-motion. In region (I) R has no real zeros. In the regions 
(II) and (V) R has two zeros. Four zeros are possible in the regions (III) and (IV). In the grey area the ^-equation does 
not allow geodesic motion. In region (b) the geodesies cross i? = but f — cannot be crossed. However, in region (d) 
r — can be crossed but $ — w is not crossed. 



Additionally one may define an effective potential consisting of the two parts V + and V by the 
equation 

R(r) = (r 2 + a 2 ) 2 (E - V + )(E - V~) , (49) 

thus 



aL± JA r [K + r 2 (5 + J 2 )] 
V ± = V , 9 - =5 -■ (50) 

The area between V + and V~ is a forbidden zone, since there R(f) < 0. V + and V~ meet at the horizons 
r±, where V ± (r±) = 1^. In the limit r ^ co and f —> — oo the effective potential V^ 1 converges to 
±v 5 + J 2 . If the sign of L changes, the effective potential is mirrored at the f-axis. 

Taking all informations into acount, we can now determine the possible orbits in the different regions 
of the parametric diagrams (below we always assume that < rj_|_i): 

1. Region (I): No real zeros and R(f) > for all r. Here only transit orbits are possible which cross 
?~ = 0. Since region (I) intersects with region (d), •& = J is not crossed. 

2. Region (II): R{r) has two real zeros r%, T2 and R(f) > for f £ (— oo,ri] and f e [ra,oo). Region 
(II) intersects with region (b) and region (d) from the i?-motion. In part (lib) there is a negative 



and a positive zero, so that escape orbits are possible for f < and two-world escape orbits are 
possible for f > 0. Here the orbits cross $ = 5 plane. The turning point of the TEO can coincide 
with the inner horizon r_. In the special case K = (aE — L) 2 the former positive zero is now at 
f = so that the two-world escape orbit turns into a terminating orbit and the motion takes place 
in the equatorial plane. 

In part (lid) R has two negative zeros, so there is an escape orbit for f < and a crossover two- world 
escape orbit. This time d = \ is not crossed. 

3. Region (III): R has four positive zeros r±, r2, r^, and R(f) > for f € [ri,T2] and f g [r3,r4]. 
Possible orbits are bound orbits with f > r + and many-world bound orbits. One or both (for L = 0) 
of the turnings points of the many-world orbit can coincide with the horizons. Since region (III) 
intersects with region (b) the orbits cross $ = 5. 

4. Region (IV): R has one negative zero r\ and three positive zero r%, r%, ri. R(f) > for f £ (—00, n], 
f G [f2, 7-3] and f € [r4, 00). Escape orbits with either f < or r > f+ and many-world bound orbits 
are possible. If K is small then also bound orbits with < f < f_ are possible instead of the 
many-world bound orbits. These orbits are hidden behind the inner horizon. In the special case 
K = (aE — L) 2 one of the turning points of the (many- world) bound orbit is now at f = so that 
the orbit turns into a terminating orbit and the motion takes place in the equatorial plane. 

5. Region (V): R has two positive zeros ri, ri and R(f) > for f G [ri,ra]. Only many- world bound 
orbits are possible. One or both (for L — 0) of the turnings points can coincide with the horizons. 
Since region (V) intersects with region (b) the orbit crosses 1? = ^. In the special case K = (aE — L) 2 
one zero is now at r = so that the orbit turns into a terminating orbit and the motion takes place 
in the equatorial plane. 

Table [2] shows all possible orbits types of orbits in the rotating black string spacetime. Some examples 
of energies corressponding to the various orbit types in the effective potential can be seen in figure [9] 

If we consider the original Kerr spacetime then the regions (III) and (V) are not present for S = 0, but 
in the rotating black string spacetime all five regions are present both for S = and 5 = 1. That means 
in the rotating black string spacetime we have bound orbits for light with f > f + (see figure 9(b) I which 
were not possible in the Kerr spacetime. 
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(a) 5 = 1, o = 0.45, L = 0.2, J = 1 and (b) 5 = 0, a = 0.4, L — 2, J — y/2 and (c) 5 = 1, a = 0.45, L = 1, J = 0.2 and 
K — 2: Examples of the orbit types A, K — 5: Examples of the orbit types D K — 0.1: Examples of the orbit type F. 
B, C and G. and E. In the case of type D, here a 

stable bound orbit for light is possible. 



Figure 9: Examples of the effective potential in the rotating black string spacetime. In grey area no motion is possible, since here 
R < 0. In the dashed area the motion is forbidden by the i9-equation. Horizontal red dashed lines represent energies and 
red points mark the turning points. The horizons are marked by vertical black dashed lines. 



3.3 Solution of the geodesic equations 

In this section we present the analytical solution of the equations of motion (l34 |) - ([38l) . 



type zeros region range of f orbit 



A Id ~ I II II TrO 

B 2 lib * I * II II EO, TEO 

B_ • — I ♦ 1 EO, TEO 

Bo * j II II EO, TO 

C 2 lid ~ ' * I II II EO, CTEO 

D 4 Bib I * II II * * * : MBO, BO 

D± 1 1 f • MBO, BO 

D_ 1 ♦ II * * MBO, BO 

D+ I * II ♦ * * MBO, BO 

E 4 TVb * I * II II * * EO, MBO, EO 

E_ •— I ♦ II * * EO, MBO, EO 

E+ * I * II ♦ * EO, MBO, EO 

F 4 TVb * I * * II II * EO, BO, EO 

F+ • — |— • •— II ♦ EO, BO, EO 

F •— + *—\\ 1 • EO, TO, EO 

Fo+ * + * \\— t EO, TO, EO 

G 2 Vb I * II II * MBO 

G± 1 1 f * MBO 

G_ 1 ♦ 1— • * MBO 

G+ 1 ♦ * MBO 

Go + i 1— • * TO 

G + + 1 f " TO 



Tabic 2: Types of orbits of light and particles in the rotating black string spacctimc. The thick lines represent the range of the 
orbits. The turning points arc shown by thick dots. The horizons arc indicated by a vertical double line. The single vertical 
line represents f — 0. 



3.3.1 The f-equation 

The right hand side of the f-equation (|34l) is a polynomial of fourth order with the coefficients: 

a 4 = E 2 -(6 + J 2 ) (51) 

a 3 = S + J 2 (52) 

a 2 = 2aE(aE — L) — a 2 (5 + J 2 ) - K (53) 

ai = K (54) 

a = & 2 [{&E- L) 2 - K] (55) 



The polynomial R = X] 4 =i can be reduced to cubic order by the substution f = ±— + tr (where ¥r is 



a zero of R): R' = Ylt=o ^ iX ' ' ■ ^ further substitution x = (Ay — transforms R' into the Weierstrafi 



form so that equation (|34[) turns into 

'' =^-9ly-gl = Pl{y), (56) 

where 

r _ H 6163 r _ hbzh b ob 2 3 b\ 

92 ~i2 1~' 93 ~^~ Te" 216- [bn 

The differential equation (|55|) is of elliptic type and is solved by the Weierstrafi p-function [57] 

2/(7) = P ~ ~fLi 92, &) » (58) 



where 7/ = 7i n + f°° > dy g with y- m = ±-£?- + If. Then the solution of acquires the form 

f = ± n — ?\ b2 + fR ■ ( 59 ) 

3.3.2 The tf-equation 

To solve the ^-equation (|3"5j) we substitute v — cos 2 $ (with v g [0, 1]) in equation (|35l) : 

d \ 2 

— ] = Av{l-v)[K -{5+J 2 )hv]-Av[hE/l-v)-L] 2 (60) 
d 7 / 

The right hand side of equation (|6T))) is a polynomial of third order 53i=i with the coefficients: 

c 3 = -Ad 2 [E 2 — (S + J 2 )] (61) 
c 2 = 4[25S(5£;-Z) -d 2 (5 + J 2 ) - K] (62) 
ci = A[K- (&E-K) 2 ] . (63) 

So equation (16"0"|) can be transformed into the Weierstrafi form using the substitution v = ^ (4w — 

^'^'-fifr = (64) 
where 

• y2 12 4 ' y3 48 216 V ' 

Equation (l6"4|) is solved by the Weierstrafi p-function and the solution $(7) of equation (|35|) can then be 
obtained by resubstitution: 



0( 7 ) = arccos (±yj ± (4p( 7 - 7 £; <?f, <?f) - |)J (66) 



With 7^ - 7in + /~ V4u 3_ d ; |B _ g g ^ = f COS 2 tf in + ff . 

3.3.3 The (^-equation 

Using the f-equation (f3~4"|) and the ^-equation (|3l))) . we can write the (^-equation (j3"6")l in the following way: 

dip = -^[(f 2 + a 2 )£ - aZl-^L 1 (aEsm 2 ^- L)^L . (67) 

A,. WR sin 2 zT y V6 

So the (^-equation consists of a f-dependent integral and a "^-dependent integral: 
V - Vin = -j-[{r' 2 + a 2 )E - aL}- 1 



I lr {(r' 2 + ~a 2 )E-~aL}^=- —^—(aEsm. 1? - L)—= = I? - I# . 

Let us first consider If. Here we substitute f — ± — ^p+rR and apply a partial fraction decomposition, 
so that If turns into 

where p\ — iSfr^j^ an d P2 — ^yzI^-Ir) 3 are ^ rst or der poles of If. The sign of ±63 depends on 
the chosen sign in the substitution r — ± — 22^— + f#. C7j are constants that arise from the partial fraction 
decomposition and depend on the parameters of the metric and the test particle. Now we substitute 
V = P (7 - 7i' n ; 52> 55) == P?(v) with v = 7 - 7i ' n : 

(70) 



The integral 1$ can be transformed in the same way by substituting first v — cos 2 then v = 
h ( 4w ~ f) and finally u = P(7 ~ Tin! at , at ) =■ Pv{v) with v = 7 - 7^: 

h = f V ~aE+ — J, dw' , (71) 

where <7 = t + f§- 

The integrals If and are of elliptic type and can be solved in terms of the elliptic p-, a- and 
£- function as shown in [51 [7] • Then the final solution of the (^-equation (|55|) is 

I \ n t \ j_ST C i ( if t \t \ , i ^f(w-«i) , CT^W + W;) 
93(7) = C (u - Win) + > , , : 2Qf{v i ){v - Win) + In — ? r - In 



1 Pf( v i) V ' Cf («in - «i) CT f (w in + Wi) 

~„ f ~ - x c 3 L ( a § {v-i q ) a#(v + v q ) \ 
ah{v - ^n) - , , /- n 2c>(w 9 )(w - Win) + m — r - in — — — — + <Pin 

4 PtfK) V ^(Win-W, ^(Win+W, / 



(72) 



where ft = pr(wi), q = ^(w g ), w = 7 - 7^, v = 7 - 7^ and 

Pf(u) =p(w,ff2>55), =P(£, 32,33) > 

0=(«) = CK <?!,<?!), =C(v,gt,gi), (73) 

cr f (w) = cr(w, 5^,55) , cr^(w) = flrf , gf) . 

3.3.4 The iw-equation 

Using the f-equation (|34[) and the ^-equation ([35]), we can write the wj-equation (|37[) in the following way: 



df di9 
dw = Jp 2 A 1 = J? 2 — + Ja 2 cos 2 ^= . (74) 

Vr Ve 



Like the (^-equation, the wj-equation consist of a f-dependent part and a ^-dependent part. We integrate 
(|74[) and use the same substitutions as in the previous section f3. 3. 31 



W - W3 in = J / Kb± 



1 , &§ 1 



2 p(w')-p I6(p(w')-P) s 



dw' + J~a 2 f (-p*(v) - dw' , (75) 
Jv in V c 3 3c 3 / 



where p = ||. The occurring elliptic intregrals of the third kind can be solved as shown in [6j [7J. 
Furthermore we use the relation J p(w)dw = — £(w). The final solution of the w3-equation is 

- ^ 1 ' -4^, - J ^f ( ^U f 2C f (w p )(w - Win) + In ^-"p) - In 

2p^(w p ) 16(p^(w p )) 3 A ^(wi„-w p ) a f (w in +w p ) 



2pf[v p ){v - Win) + 2(Cf (w) - Cr(Vm)) + 



!6 (pf^p)) 2 V "" Pf ( u ) - Pf(v P ) pf (v in ) - pf(v p ) 

(Q(V) ~ Q(V- m )) 5 (W - Win) + Wi n , 

C3 3c 3 

(76) 

where p = p?(v p ). 
3.3.5 The t-equation 

Using the f-equation ([3^)1 and the ^-equation ([33)) , we can write the f-equation (|38p in the following way: 

d < = r l±^L [(f 2 + a 2 )£; - fiZ] 4^ " «(^ sin 2 j? - L) • (77) 

Like the (/9-equation and the w-equation, the ^-equation consist of a f-dependent part and a ^-dependent 
part. We integrate (1771) and use the same substitutions as in section [5.3.31 

t-Un = f [c' + ^ — ^ + ^ ) dw' - f fa 2 i? - aL + % - —p*{v)) df , 

(78) 



where pi = ^j^fe—^ff 3 an d Pi — ^^(f ^-r) 3 are nrst 0T( ^ eT P°l es an d P3 = y§ is a second order pole. 
C[ are constants arising from a partical fraction decomposition of the f -dependent part. The occurring 
elliptic intregrals of the third kind can be solved as shown in [5J |7]- Furthermore we use the relation 
/ p(v)dv = — £(v). The final solution of the f-equation is 

n \ nit \ , sr^ C'i ( n + { \ i n , , cr?(v-Vi) a f (v + v{) 



■(Pf( v i)\ " CXr(v m — Vi) (Jf(v in + Vi) 

C 3 7-77^ \ \ 2 I 2Cf (w 3 )(« - ^n) + In — — - ln 



(Pf( w 3)) 2 V " (Tr{Vin-V3) CTr(v in + V 3 ) J (79) 

2pf (w 3 )(w - U in ) + 2(Cf (Uj - Cf ('Win)) + 



(pfM) 2 V "' Pf(v)-Pf(«3) Pf(Uin) - Pf («3) 

/ 2 C2\ 4a 2 i? 

- ( a £/ - aL + — (« - Win) (Cl?W ~ ( u in)) +^in , 

V 3 / C3 

where p 4 = pf{v t ). 



3.4 The orbits 

With these analytical results we have found the complete set of orbits for light and test particles in the 
rotating black string spacetime. Depending on the parameters 5, L, J, a and E various orbits are possible. 
Note that the horizons in the Kerr spacetime and in the rotating black string spacetime (x-y-z-plots) are 
not spheres but ellipsoids. A x-y-w-^lot of the horizons in the rotating black string spacetime will show 
cylinders. 

Figure [10] shows a terminating orbit in the rotating black string spacetime, in the x-y-z-plot this orbit 
lies entirely in the equatorial plane. To hit the ring singularity at p 2 = the parameters corresponding to 
the orbit have to fullfill the condition K = (aE — L) 2 . Otherwise it is possible for a geodesic to cross f = 
without touching the singularity. Like in the Kerr spacetime the radial coordinate f can take negative 
values. A transit orbit wich starts at positve f, crosses f — and the continues at negative r can be seen 
in figure [TT] When f changes from positive to negative it looks like the particle is reflected, which can be 
interpreted as gravity becoming repulsive for negative f- values (see [251 1251 [5D] V 

An escape orbit is shown in figure 1121 If the turning point of an escape orbit lies behind the two 
horizons the orbit is called two- world escape orbit (see figure ITS")) . Since both horizons are traversed twice, 
the test particle or light emerges into another universe. It is also possible that the turning point is at 
negative f (see figure HH) , then the orbit is called crossover two- world escape orbit since it does not only 
traverses both horizons, but also crosses f = 0. 

A special feature of the the rotating black string spacetime are stable bound orbits of light (see figure 
I15[) . Such orbits do not exist in the ordinary four-dimensional Kerr spacetime. 

Like in the Kerr spacetime there are also bound orbits hidden behind the inner horizon f_ if the Carter 
constant K is rather small (see figure [T6|) . 

It is also possible that a particle or light on a bound orbit crosses both horizons several times. Each 
time both horizons are traversed twice, the orbit continues in another universe. Such a many- world bound 
orbit is depicted in figure ITT1 

Figure [T51 shows a so called "zoom- whirl" bound orbit (see also 'A 1 32, 33 ). Here the energy of the test 
particle is very close to a local maximum of the effective potential, so that the particle "whirls" around 
the black string before reaching the periastron and then "zooms" out to return to its elliptical orbit where 
it moves to the apastron. 



(a) x-y-z-plot 



(b) x-y-w-p\ot 



Figure 10: S = 1, a = 0.45, L = 0.8, J = 0.2, K = (aE - L) 2 = 2.1025 and E = 5: 

Terminating orbit for particles in the rotating black string spacctimc. The ellipsoids or cylinders arc the horizons. In the 
left picture, the position of the ring singularity is marked by a red circle. 




(a) x-y-z-plot (b) x-y-w-p\ot 



Figure 11: 5 = I, a = 0.4, L = 0.6, J = 1.5, K = 1 and E = 5: 

Transit orbit for particles in the rotating black string spacctimc. The ellipsoids or cylinders arc the horizons. In the left 
picture, the position of the ring singularity is marked by a red circle. 



(a) x-y-z-p\ot (b) x-y-w-p\ot 

Figure 12: 5 = I, a = 0.45, L = 1.5, J = 0.6, K = 2 and E = 1.251: 

Escape orbit for particles in the rotating black string spacctimc. The ellipsoids or cylinders are the horizons. 



(a) x-y-z-plot 



(b) x-y-w-p\ot 



Figure 13: 5 = 0, a = 0.45, L = -0.5, J = 2, K = 5 and £7 = 2.25: 

Two-world escape orbit for light in the rotating black string spacctimc. The ellipsoids or cylinders arc the horizons. 





ft: 



(a) x-y-z-plot 



(b) x-y-w-p\ot 



Figure 14: 5 = 1, a = 0.45, L = 0.8, J — 2, K — 2 and £ = 5.5: 

Crossover two-world escape orbit for particles in the rotating black string spacctimc. The ellipsoids or cylinders arc the 
horizons. 





(a) a:-y-2;-plot 



(b) x-y-w-p\ot 



Figure 15: 5 = 0, a = 0.4, L — 2, J — >/2, K = 5 and E = 1.36: 

Bound orbit for light in the rotating black string spacctimc. The ellipsoids or cylinders arc the horizons. 



(a) :c-y-.2-plot 



(b) x-y-w-p\ot 



Figure 16: 5 = 1, a = 0.49, L = 1.4, J = 0.5, K = 0.6 and E = 1.63: 

Bound orbit for particles hidden behind the inner horizon in the rotating black string spacctimc. The ellipsoids or 
cylinders are the horizons. 





(a) x-y-z-plot (b) x-y-w-plot 

Figure 18: 5 = 1, a = OA, L — 2, J — V% K = 5 and E = 1.706: 

"Zoom-whirl" bound orbit for particles in the rotating black string spacctimc. The ellipsoids or cylinders arc the horizons. 



4 Conclusion 



In this article we presented the complete set of analytical solutions of the geodesic equations of test 
particles and light in the static and the rotating black string spacetime. The static and the rotating black 
string metric are obtained by adding a compact dimension to the Schwarzschild and Kerr metric. 

The geodesic equations can be solved in terms of the elliptic p-, a- and ^-functions. Using effective 
potential techniques and parametric diagrams, possible types of orbits were derived. In the static case 
BOs, EOs and TOs are possible, while in the rotating case BOs, MBOs, EOs, TEOs, CTEOs, TrOs and 
TOs are possible. 

In contrast to the ordinary four-dimensional Schwarzschild and Kerr spacetime, bound orbits of light 
are possible both in the static and the rotating black string spacetime. 

Hereby the analytic solutions do not only give a proof of the existence of bound orbits of light in the 
black string spacetime. They also present a usefull tool to calculate the exact orbits and their properties, 
including observables like the periastron shift of bound orbits, the light deflection of escape orbits, the 
deflecton angle and the Lense-Thirring effect. For the calculation of the observables analogous formulas 
to those given in [12] can be used. Observables can later be compared to observations. 

The black strings considered here are so-called uniform black strings, since there is no dependence 
of the extra dimension. However, uniform black strings are subject to the Gregory-Laflamme instability 
[34) . which is associated with the emergence of a branch of non- uniform black strings, whose horizon size 
is not constant w.r.t. the compact direction, but depends on the compact coordinate [35l 023 [37]. At 
the endpoint of this nonuniform black string branch a horizon topology changing transition should be 
encountered [551 I3"9"] . 

To gain a better understanding of this transition, it is essential to solve the geodesic equations. How- 
ever, for the construction of nonuniform black strings no analytic techniques are available, and numerical 
techniques have to be employed. Likewise, the geodesic equations must be solved numerically. It is there- 
fore of high relevance, to have analytic solutions of the geodesic equations available as a testbed. Moreover, 
it will be very interesting to see, how the set of analytic solutions changes, as the black strings become 
more and more nonuniform. 

Also non-uniform rotating and charged black strings have been considered in |40j and |41j whose 
geodesies would be very interesting to study. 

It is also interesting to consider higher dimensional (d > 4) Schwarzschild spacetimes, here stable 
bound orbits are no longer possible [11] . One may wonder wether bound orbits of light become possible 
if a compact dimension is added to the metric. But it can be shown that unlike in the (rotating) black 
string spacetime in five dimensions, stable bound orbits of light are not possible if a compact dimension 
is added to higher dimensional Schwarzschild spacetimes. 

Moreover the bound orbits of light appear neither in the spacetime of an Abelian-Higgs string [23] nor 
in the spacetime of two interacting Abelian-Higgs strings [24] . There escape orbits are the only possibility 
for massless particles. For cosmic superstrings bound orbits of light are possible in Melvin spacetimes but 
not in asymptotically conical spacetimes [25) . 
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